Evolution of Parton Fragmentation Functions at Finite Temperature 



Jonathan A. Osborne^, Enke Wang^'^ and Xin-Nian Wang^'^ 
^Nuclear Science Division, Lawrence Berkeley Laboratory, Berkeley, California 94720, USA 
^Institute of Particle Physics, Huazhong Normal University, Wuhan 430079, China 
^Department of Physics, Shandong University, Jinan 250100, China 
(October 25, 2002) 

The first order correction to the parton fragmentation functions in a thermal medium is derived 
in the leading logarithmic approximation in the framework of thermal field theory. The medium- 
. . . modified evolution equations of the parton fragmentation functions are also derived. It is shown 

OA , that all infrared divergences, both linear and logarithmic, in the real processes are canceled among 

■ themselves and by corresponding virtual corrections. The evolution of the quark number and the 

O . energy loss (or gain) induced by the thermal medium are investigated. 
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I. INTRODUCTION 



Parton production in hard processes is normally followed by final-state radiation and subsequent hadronization, 
. giving rise to collimated jets of hadrons as observed in experiments. The hadron distributions inside a jet known as 
' jet fragmentation functions can be defined as the vacuum expectation values of the parton fields and the hadronic 
. interpolating fields Though these fragmentation functions are non-perturbative and currently can only be mea- 
I sured experimentally, their evolution with the probing scale can be calculated within perturbative QCD. The resultant 
Dokshitzer-Gribov-Lipatov-Altarelli-Parisi (DGLAP) |^ evolution equations have been stringently tested against ex- 
' periments and now can even be used to measure the scale-dependence of the running strong coupling constant |3j . If 
[ the hard parton is produced amid a thermal QCD medium, as most likely occurs in high-energy heavy-ion collisions, 
f-^ ^ the subsequent final-state radiation and parton cascade must be modified by the presence of the medium. One not 
I only has to consider parton emission but also parton absorption Q in describing the evolution of the fragmentation 
functions. Such a medium effect is closely related to the radiative parton energy loss, which also leads to modified 
fragmentation functions In most of the studies of parton energy loss, one relies on a Debye-screened potential 
model [^,0j^J^ for parton-mcdium interaction. One then has to introduce additional parameters such as the mean free 
path, which is ill-defined when the magnetic part of the one-gluon exchange interaction is included. One can replace 
the mean free path by a transport parameter Q and eliminate the infrared problem by using a hard-thermal-loop 
; ^ (HTL) resummed gluon propagator ]lO[ |. However, it is still theoretically interesting to study the problem within the 
framework of field theory at finite temperature. 

This paper is our first attempt to study medium-modified fragmentation functions within QCD field theory at 
finite temperature. We will first extend the definition of the fragmentation functions to include the scenario of parton 
propagation inside a thermal bath. We will calculate the first order corrections to the parton fragmentation functions 
in a thermal medium in the leading logarithmic approximation and derive the corresponding modified evolution 
equations. We will find that the modified evolution equation can be cast into a similar form as the DGLAP equations 
in the vacuum. However, the modified splitting functions depend explicitly on the temperature and the partons' initial 
energy. We will study the structure of different contributions and their physical interpretations. We will demonstrate 
that all infrared (both linear and logarithmic) divergences in radiative corrections cancel either among themselves or 
with the virtual corrections. We will also study the evolution of the net quark number and energy loss (or gain). 

We should emphasize that our calculation in this paper includes only the first order corrections in the leading 
logarithmic approximation. At this order, we can only consider parton emission and absorption. By solving the 
evolution equations, one can effectively resum radiative corrections associated with the leading parton. However, 
collision-induced radiation is not considered at the first order of perturbation. Furthermore, we cannot include any 
interference effect like the Landau-Pomeranchuck-Migdal (LPM) interference pd[ | in medium-induced bremsstrahlung 
in the leading log approximation. One has to include the spectral function of the HTL resummed gluon propagator 
and go beyond the leading log approximationin order to consider radiation induced by multiple scattering and the 
LPM interference effect. We leave this to future investigation. 

The remainder of the paper is organized as follows. In Sec. II, we will first review the definitions of the fragmen- 
tation functions in the vacuum and the basic physical processes that lead to the DGLAP evolution equations of the 
fragmentation functions. We then extend the definition to include the case in which a parton propagates through a 
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thermal QCD medium. In Sec. Ill, we derive the radiative corrections to the fragmentation functions in a thermal 
medium based on finite temperature field theory. We will discuss various physical processes that contribute to the 
radiative corrections and the cancellation of all infrared divergences. In Sec. IV, the QCD evolution equations are 
derived from the calculated radiative corrections. As examples of the application of the evolution equations, we derive 
the evolution equations for the net quark number and calculate the energy loss (gain) induced by the medium. Finally, 
we give a summary and conclusion in Sec. V. 



II. FRAGMENTATION FUNCTIONS AT ZERO AND FINITE TEMPERATURE 



A. Fragmentation Function at Zero Temperature 

Semi-inclusive cross sections can generally be factorized into the convolution of a parton fragmentation function 
and a hard partonic cross section in the coUinear (leading twist) approximation. Fragmentation functions are defined 
as the vacuum expectation values of parton fields and hadronic interpolating fields. The e'^e~ annihilation process 
is an ideal framework in which to define the quark and anti-quark fragmentation functions. The total hadronic cross 
section for this process can be expressed as 

CTe+e-^X = —^L^,{pa,Pb)W''-'{q) , (1) 

s 

where s = (pa + Pb)'^ is the invariant mass of electron-positron system and q = Pa + Pb the momentum of virtual 
photon (7*). The leptonic tensor L^^ is given by 

L^u{Pa,Pb) = ^Tr(7^^„7^^b) , (2) 
where 1/4 comes from the average over spin polarization of the initial e+e^ state. The hadronic tensor is defined as 

^"^('^) ^^T.(o\j''iom{x\riomi27rr6\p,-q) 

X 

= ^ y"d'*xe-*«-^(0|J^(0)J''(x)|0) , (3) 

where ^'^ns over all possible intermediate states and the quark electromagnetic current is = '^gQq4'qlfj.'4'q- 
Here, Qq is the electric charge of the quark in imits of the proton charge. 

At the lowest order in pQCD, the electron-positron pair annihilates, forming a virtual photon which then decays 
into a quark anti-quark pair. After a sufficiently long time, the partons undergo a nonperturbative fragmentation 
process and emerge from the scattering center in the form of hadronic jets. The two partons fragment independently 
of one another in the leading twist approximation, allowing us to define process-independent fragmentation functions. 
In terms of these fragmentation functions, the semi-inclusive cross section Ue+e-^h can be shown to have the form 

= cyf [D,^u(zh) + Dq^nizn)] (4) 
azh ^ 

in the center-of-mass frame. Here, 

of ^ Nc^^ (5) 

is the perturbative cross section in the lowest order for e^e~ qq and D q(^q-^^h{zh) is the probability that a quark 
(anti-quark) will decay into a hadron h with a fraction Zh of its energy. is the number of quark colors. The energy 
of the quark (anti-quark) is just ^'^(g) — \fsl'2., so the observed hadron has energy — Zh\/sl2. This relationship 
can also be written in the invariant form 

Zh = 5 — • (6) 

q 

Applying the coUinear approximation to Eq. dST), one can obtain a formal definition of the fragmentation functions 
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where the Dirac operators Tq^q){k,ph) are given by 

{fq)o.p{k,Ph) = j d^xe-'''-^Y.^Q\',i;^{Q)\ph,S){ph.S\i:p{xm 



(7) 



(T^)aMfc-P'«) = / d*xe-*-^(0|V;/5(0)K,5)(p,„^|^„(a;)|0). (8b) 



Here, the sums go over all physical states S and sums over quark colors are implicitly averaged. The light-like vector 
n*^ = [n^ , n~ , n±\ = [0, 1, Ox] is taken conjugate to the momentum of the observed hadron in the sense that its spatial 
components are antiparallel to the spatial momentum of the hadron. This implies n ■ ph = p^ = {p\ + |p/i|)/2. The 
gauge links required to make this expression gauge-invariant have been suppressed since they do not contribute to 
the leading-twist fragmentation functions in the light-cone gauge, n ■ A — Q. 
The gluon fragmentation function Dg^fi{zfi) are defined in a similar way: 
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5(zu~'^-^)d^Ak)f^'^{k,ph). (9) 



Tr{k,Ph) = / d4xe-^'=-^^(0|A''(0)K,5)(p;„5|^^(x)|0). (10) 
•' s 

Similarly, gluonic color indices are implicitly averaged over. The gluon polarization vectors £^{k) satisfy 

dt,y{k)= ^ ef,{k,X)£i,{k,\) 

A=l,2 

kfiTiu -f- kiiTifi n 1 A 

" ^^""^ + ~k ^^^> 

in the light-cone gauge. The field-strength tensor, G^i, — d^Ay — d^A^ + ig[A^^ A^], allows us to write the gluon 
fragmentation function (^) in the manifestly gauge-invariant form 

D,->h{zn)^^^Y. J ^e-^^"/^''(0|G+'^(0)K,5)(p.,5|G%(x-)|0). (12) 



s 



Once again, the gauge links necessary to make this expression gauge-invariant do not contribute to our functions in 
the light-cone gauge. 

Using the 'cut vertices' technique introduced by A. Mueller |l^, we can calculate the scale dependence of these 
functions in perturbation theory. The Feynman diagrams illustrating this calculation are shown at leading order in 
Fig. The corresponding bare quark and gluon cut vertices are 

o ^r'' — z — j 

Zph ■ n \ k ■ n / 

and 

Zhdf,t,{k)S (^Zh - , (13b) 

respectively, in the light-cone gauge. 

At next-to-leading order, 0{as), the gluon bremsstrahlung process represented in Fig. |^ gives a positive contribution 
to the quark fragmentation function. This is due to the fact that more options are open to the fragmenting quark 
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at this order: it can fragment immediately, or radiate a gluon first. This bremsstrahlung process becomes easier 
as the energy of the radiated gluon decreases, leading to a divergence as the energy goes to zero. Physically, this 
divergence is meaningless. If the quark radiates no energy, its fragmentation should not be altered. This idea is 
realized mathematically by the self-energy contribution represented in Fig. ^j. This correction removes the double- 
counting inherent in the above argument, and ensures that each decay channel is counted only once. For this reason, 
the divergent contribution from zero-energy bremsstrahlung is entirely canceled. 

All next-to-leading order fragmentation sub-processes lead to similar corrections to the fragmentation functions. 
We will list a complete set of Feynman 'cut' diagrams later. Apart from the spurious infrared divergences in some 
of the diagrams, each correction contains a coUinear divergence as the transverse momentum of the loop approaches 
zero. These divergences signal an unavoidable scale dependence in our corrections. Physically, this scale dependence 
represents the fact that partons appear different when probed at different scales. We can handle these divergences 
mathematically by computing only the difference between our fragmentation functions probed at different scales. Since 
the asymptotic contribution to our loop integrals cannot depend on the external scales of the process, the collinear 
divergences will always cancel in such differences. With this in mind, we write the results 
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of the calculation. Here, we have defined the singlet quark fragmentation function as. 



,h{z,li^) = ^ {Dg^h{z,^i^) + Dg^h{z,fl^)) 
1 



(15) 



The integration kernels, or splitting functions, are given by [|l3| 
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Here, is the number of active quark flavors and the SU(iVc) Casimirs are given by Cp = {N^ 
and Tp = 1/2; = 3. The '-|- '-function is deflned such that the replacement 
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is valid for any function f{z) that is continuous at z — 1. 

The scale dependence of the fragmentation functions is then determined by the DGLAP evolution equation 
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^qq{z)Dq^ll{zh/z, Q^) + Jqg{z) D g^h{zh / Z , Q^) 



Jgq{z)Ds^h{Zh/z, Q^) + "fgg{z)Dg^h{Zh/ Z, Q^) 



(18a) 
(18b) 



as can be shown by repeated iteration of two-particle-irreducible diagrams similar to those shown in Figs. ^ and 
Eq.(p^ is analogous to the 'Born' approximation of (|l^). Since we have calculated our splitting functions only at 
leading order in a^, our evolution equations are not correct to all orders. However, it can be shown |l|] that Eq. (|l|) 
contains the highest power of log(Q^//x^) present at each order in as- For this reason, it is known as the leading log 
approximation (LLA) to the full evolution equation. 
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of the zero temperature functions cannot be complete. In particular, fragmentation functions relating to different 
initial parton energies will mix with each other under renormalization, resulting in a very complicated evolution 
equation. 

A more straightforward way to approach the problem of fragmentation at finite temperature is to consider the 
momenta of the fragmenting parton and the final state hadron separately. In particular, 

Fk/,{y;x) = ^ I -^S{k+^x) I d^ze-"'^J2T^T[{m\S,PH)l+{S,p,Mz))] (22) 

represents the probability for a quark of momentum p"*" = x to decay into a hadron of momentum p'^ = y. This 
expression is to be understood as a thermal average over the quark field operators in which a particular intermediate 
state, \ph), is singled out. The set of states {\S)} includes the thermal phase space indicated in Eq. ( |l9| ) and the 
quark color is implicitly averaged. Note that in the following discussion of the newly defined fragmentation functions 
at finite temperature, all variables, x, y, z, etc. represent absolute '+'-momentum and thus carry the dimension of 
energy. We will denote fractional momentum by these variables with a subscript {zf, for example). 

This new function is quite similar in form to the quark fragmentation function in the vacuum. Defining 

Dq^h{zf;x) = F,,/q{xzf;x) , (23) 

we find that Dq^h{zf; x)\t=q — Dq^hizf). However, the interpretation oi Ff^/q is entirely different. By tracking the 
absolute momentum of both the initial parton and the observed hadron, this new function is able to clearly express the 
dependence of fragmentation on the temperature and the initial parton energy. We will see below that this separation 
of fragmenting parton and observed hadron is essential to a transparent evolution equation at finite temperature 
where fragmentation functions with different initial parton energies mix. 

Fragmentation functions for anti-quarks and giuons can be defined in a similar way. We have 

Fn/Tiiv^ ^) = f / 5{k+ -x) ( d^ze-^'^' ^(V^(0)|5,p,)7+(^,P.|V(^)) (24a) 
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Fk/g{y;x) = -^ jSlSik+~x) /'dSe-'^-^(G+''(0)|^,p.)(^,p.|G+(z)) , (24b) 

s 

where colors are implicitly averaged over, as before. 

III. RADIATIVE CORRECTIONS 

In this section we calculate the first order correction to the fragmentation functions at finite temperature. The 
result will be used to derive the corresponding evolution equations in the next section. Using 



ySik+~x) = zl^s(z,-Pf^) , (25) 
X ph ■ n \ k ■ n J 



where Zh = y/x and ph-n ~ y, we see that the form of this calculation is quite similar to that of the zero-temperature 
case. The only practical difference between the two cases lies in the form of the propagators, both cut and uncut. For 
uncut (intermediate state) propagators, the expressions 

^(^) = - 2^/(fc+)^(fc')l , (26a) 

^M-(fc) = \t2^ + 2™(fc+)5(fc2)l d^,{k) , (26b) 
replace the ordinary zero-temperature propagators. The distribution functions 

/(^) - ' ^^^^^ 
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Using the collinear approximation, 



l = h=[vlz'f,Q,0^], 

in Eq. decouples the integral over I in Eq. (30) from H{lc)- This allows us to write the correction as 



(32) 



xz, 



Performing the change 



in H(lc), we arrive at the result 
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Here, we have defined 
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dz 
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The integral over z in Eq.(|35D can be divided into two regions: < z < a; and z > x. In the first region, we have 
the physical situation in which the quark emits a gluon into the thermal bath before it fragments. The second region 
represents a physical absorption process. In either case, the plus-component of the intermediate quark momentum is 
z. At zero temperature, only the emission process contributes. Since the quark can only lose energy through evolution 
in this case, the '+ '-momentum of the intermediate quark must be larger than that of the observed hadron as well as 
smaller than that of the initial quark, y < z < x. This leads to the limits on the zero-temperature evolution equation 
in Eq. (18). At finite temperature, the quark can gain energy through absorption. This process extends the region 
of integration in both directions, z > x and z < y. The extension to z > x is an obvious effect of absorption: the 
intermediate quark in an absorptive process has more energy than the initial quark. The extension z < ?; is somewhat 
more subtle. At zero temperature, conservation of energy implies that the observed hadron cannot have more energy 
than its parent parton. At finite temperature, we have no such restriction. The thermal bath can provide the excess 
energy needed for a parton to fragment into a more energetic hadron. This can be thought of as a kind of 'collective' 
fragmentation of the parton and the bath. We can immediately see from Eq. (^5|) that the inclusion of this effect 
is necessary: even if we choose to take Fh/q{y\x) = for y > x, Eq. ( ^sf ) will generate a finite contribution in this 
region. 

Our phenomenological interpretation of Eq. ( p5[ ) leads to a somewhat fundamental difficulty with the definition 
of fragmentation functions at finite temperature. As the intermediate energy, z, decreases around and below the 
temperature, it becomes inconsistent to treat the parton as separate from the bath itself. Furthermore, it is impossible 
to distinguish the fragmentation of a low-energy parton in a thermal bath from the fragmentation of the bath itself. 
This difficulty can be resolved via a phenomenological separation between the fragmentation of a parton distinct from 
the bath and that of the bath itself, as we will discuss in the next section. For the remainder of this section, we will 
ignore our difficulty and take the radiative corrections at face value. 

As in the zero-temperature case, our result for real emission/absorption is plagued with infrared divergences as 
z —> X. The divergences in Eq. ( |35[) are made worse by the Bose-Einstein functions. As z ^ x, these functions 
generate a linear divergence in addition to the logarithmic divergence present in the zero-temperature part. All of 
these divergences are unphysical in the sense that they represent overcounting of the relevant degrees of freedom. As 
explained above, they must be canceled by virtual corrections. 

The virtual correction shown in Fig. reads 
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Here, £(/) is the quark self-energy. The prefactor of 1/2 comes from the renormahzation of the initial state ||l6|| . 

As before, we intend to use the coUinear approximation in Eq. ( ^ ) to relate the nonperturbative part of this 
expression to i^/i/q- In this case, however, direct substitution of I — Ic leads to an indeterminate form, since Tt{l)/P 
is not defined for P — l^. For this reason, we must carefully take the limit of Eq. ( p9| ) as / Ic- With this in mind, 
we can read the quark self-energy at one-loop from the diagram in Fig. ^ 
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After completing the k integral by contour, the coUinear divergent part of this expression is given by 
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dz 
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°° dz 



Pgg{z/x) [(1 + n{x ~ z)) f(z) + n{x - z) (1 - f{z))] 



Pgg{x/z) [n{z -x){l- f{z)) + (1 + n{z ~ x)) f{z)] 



(41) 



where we have taken I = Ic wherever it does not generate a singularity and used the property Pqq{zf) = —ZfPqq{l/zf) 
of the quark splitting function. This result can be understood in terms of physical decay and scattering processes at 
finite temperature |1^. In the second equality above, the first line corresponds to the decay of a quark into a quark 
and a gluon, as well as the inverse process. The second and third lines correspond to the forward and backward 
processes of Landau damping via scattering of the quark with another thermal quark or gluon from the thermal bath. 
Landau damping requires the presence of thermally excited states, and disappears at zero temperature. In each term, 
every incoming gluon (quark) from the thermal bath is weighted with a Bose-Einstein distribution function n{kf /T) 
[a Fermi-Dirac distribution function fikf/T)], and each outgoing gluon (quark) receives a Bose-Einstein enhancement 
factor 1 + nikl/T) |aPauh blocking factor 1 - f{kf /T)]. 

Substituting Eqs.(p9|) and ( ^ ) into (^) and collecting terms, we arrive at the virtual correction 



AF(;)(2;;a:) = -gF,/,(y;x) j ^ p^'^[e{z)e{x - z) + e{x - z)n{x - z) - e{z)f{z)]Pqq{z/x) . 



(42) 



In order to see clearly that all infrared divergences in the real gluon radiation are canceled by the virtual correction, 
we combine Eqs.(|3^) and (E2|) and divide the integral into different regions: 
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dz z 
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As z X, the logarithmic divergences in terms (a), (6) and (d) are canceled by those in (c), (/) and (g), respectively. 
The linear infrared divergences in terms (a) and (b) are canceled by those in (c) and (/), leaving behind a residual 
logarithmic divergence. This divergence cancels in the sum of terms (a) and (6) as long as is differentiable 

at z = X. Term (e) is finite. 

Another real correction to quark fragmentation comes from Fig. ^, 



X / d4^e-*'«^(G+'^(0)|^,p,)(5,p,|G+(e)) d.^mi^'^^^Hl) : 



(44) 



Here, the loop gives the hard contribution 
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Performing the integrals and simplifying, we obtain 

dfc^ f°° dz 



The new splitting function, 
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(45) 
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(47) 



is related to Pqq(zf) via z/ — > 1 — z/. 

As in the zero temperature case, the scale dependence of the fragmentation functions is generated by the coUinear 
divergences present in Eqs. (^ ) and (^6|). To remove the divergences, we consider the difference between fragmentation 
functions measured at two different scales: 



Fh/g{y; X, Q-') - F^/qiy; x, //') = ^Fjf/^iy, x, Q'') + ^Fl)[{y; x, Q'') + ^Fl^,[{y- x, Q^) 



(9), 
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lqq{z;x)Fh/q{y]z,y?) +-fqg{z;x)Fh/g{y]z,^i^)\ . (48) 



The modified kernels are 



^qq{z\x) = Pqq{z/x)e{x - z)[9{x - z) + n{x - z)] 

/OO 7 
— Pqq{io/x) [e{uj)e{x - Lu) + t{x - uS)n(x - uS) - e(cj)/(w)] (49a) 
-OO X 

^qg[z;x)^Pqg{zlx)t(x- z){Q(x- z)- j{x~ z)\ , (49b) 

where the splitting functions are given by Eqs.(|3^) and (p7|). Obviously, the zero-temperature result is recovered as 
n and / approach zero. 
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Simplifying, we obtain 



d^p{l)Vf'^''{l, k-l, -k)^^^-!, I - k, k)dxAk - l)df,p{k)dP^ {k)dr,{k) 
= d^p{l)V^^''{l, k^l, -k)V"^-'{-l, I - k, k)dxa{k - l)d^,{k) 



where Zj = y/xz'^ as before. Performing the k^ and k integrations leaves us with 

AF(^)(y;^) = g / '^P,,{z /x)e{x ~ zMx - z) + n{x - z)]F^,g{y- z) ^ 

The splitting function 
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is familiar from the zero- temperature result. 

Ks z ^ X, our correction in Eq. ( |6l| ) generates infrared divergences which must be canceled by virtual contributions. 
The loop in the diagram of Fig. p| represents 
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n^'^H) is the non-Abelian contribution to the gluon self-energy, which can be written as 



nr(0 = Ig'CA 



d-^k 



Vf^^^ik, I - k, ~l)V°'^%-k, k-l, l)Ao,0(k)Ax^{l - k) 



(64) 



The symmetry factor 1/2 takes into account the indistinguishability of the intermediate gluons. As with the quark 
virtual correction, we must be careful when taking the limit I ^ Ic- The result of the calculation contains the coUinear 
divergence 



Ihn 1 d,^{l) d,p{l)dP''{l) - ^hn 1 [n^''(0] d,^{l) 



27r 



dz 



/ - - ^) + 2n(z)e(z)] P,,(z/a;) / 



dki 



kl 



(65) 



As in Eq.(p4|), Wi'^ {l)dy^j,{l) /P can be divided into physical decay and scattering processes coming from Landau 



damping at finite temperature. This leads to the contribution 

1.2 /-oo 



AF(;f(y;x) = -^F,/,(2/;x) j ^ pJ-^[e{z)e{x~z) + 2e{z)n{z)]P, 



,(z/x) 



(66) 



Dividing this integral into different regions, as in Eq. (^), one can see that all of the infrared divergences present in 
Eq. ( |6l| ) are canceled by those in Eq. (|66|). 

An additional virtual correction to the gluon's fragmentation function is illustrated in Fig. ||. Although this 
correction contains no infrared divergence, it has the normal coUinear divergence and thus will contribute to the scale 
dependence of the fragmentation function. This correction contains the hard part 



nr(07T^+nr(0- ' 



P + it] 



P — irj 



ld,Al)d^-il)5(y~^ 
x \xi+ 



n/Im 



P -\- irj 



!Ld„r{l)d^'^{l)5[^--j^ 



y pI 



(67) 
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Here, n^^(/) is the quark contribution to the gluon self-energy. It has the form 

nr (0 = 9sTf / Tr [j^S{k)rS{k (68) 



Standard manipulations lead to a coUinear divergent term 



;hn i d^^il) [n^/il)] d,rm%l) = Imi i [W^%l)d,M 

= / -[e{z)e{x -z)- 2f{z)e{z)]P,,{z/x) / , (69) 

which gives rise to the virtual correction 

AF(;^)(y;x) = -^F,/,(y;x) j ^ p^'^[e{z)e{x - z) -2e{z)f{z)]P,,{z/x) . (70) 

This correction does not contain an infrared divergence, as mentioned earlier, but the presence of a coUinear divergence 
implies that it will still contribute to the scale dependence of the fragmentation function. 
Putting all of these corrections together, we arrive at the expression 

Fh/Av^ X, Q") - Fh/,{y- X, fj,') = AFjf/^y; x, Q^)+AFl%{y; x, Q')+AF^;/J (y; x, Q')+AF^;J (y; x, Q^) 

= 27r y 2 y —['laq[z;x)Fh/s[y■,z,^J.) + Jgg(z■,x)Fh/g(y■,z,^J.)\ (71) 

for the renormalized gluon fragmentation function. The modified kernels are 

jgq{z;x) = Pgq{z/x)e{x - z)[9{x - z) - f{x - z)] (72a) 
^gg{z;x) = Pgg{z/x)e{x - z)[6{x - z) + n{x - z)] 
1 f^d 

~-5{l~xlz) / —{Pgg{^JJ/x)e{x-Lo)[e{u)e{x-^JJ) + 2n{x-Lo)] (72b) 
+2nfPgq{uj/x)e{uj)[9{uj)9{x - w) - 2/(w)]} , 
where the splitting functions Pga{zf) have been given in Eqs.(|5^) and (|62|). 



IV. EVOLUTION EQUATIONS AND ENERGY LOSS IN A THERMAL MEDIUM 

The first order corrections to the fragmentation functions derived above can be used to resum the leading powers of 
as{n'^) log((5^//i^) in our perturbative expansion. In order to obtain these leading logarithms, the transverse momenta 
of the emitted and absorbed thermal partons must be ordered from largest to smallest as the quark propagates. 
Essentially, this fact arises from the form of the nested loop integrals. In order to obtain a double logarithm from 



dfc2 fQ- dll 



(73) 



we must have > k'j^. At zero temperature, this ordering implies that the virtuality of the leading parton is similarly 
ordered. Absorption processes present in thermal media disallow this virtuality ordering, so only the minimum 
requirement of transverse-momentum ordering is observed by our corrections. 

As in Eq.(|l8|), the result of the resummation can be expressed in the compact form of an evolution equation: 

Q'^^^V9(y;2;,Q') = ^^y^ ^\lgq{z-x)F^I,{y-^M'')^l9a{z-x)F^lg{y-zM^)\ ■ (74b) 
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These equations express the change in our fragmentation functions as the probing scale is varied. At zero temperature, 
this dependence arises from the arbitrary distinction between a parton and its surrounding vacuum fluctuations. The 
theoretical process in which a quark emits a gluon with transverse momentum smaller than the probing scale cannot be 
distinguished from the process in which no such emission occurred. Hence, this state is considered part of the vacuum 
fluctuations inherent in the nature of the quark's propagation. On the other hand, if the transverse momentum 
between the quark and emitted gluon is larger than the probing scale, the final state is physically different from that 
of a single quark. In this case, the process represents a physical bremsstrahlung. Obviously, changes in the probing 
scale will give rise to changes in the effective fragmentation function. 

At finite temperature, the probing scale also serves to separate the parton from the surrounding thermal bath. Just 
as we cannot completely separate a parton from its vacuum fluctuations, a parton in a heat bath cannot be completely 
separated from the fluctuations of the bath. The idea behind these two cases is quite similar, but there are some 
striking qualitative differences. In particular, the vacuum does not carry any energy. Therefore, the total energy of a 
fragmenting parton is conserved under changes of scale: 

^2 ^ ^ / dzzD,^hiz,Q^)=0 . (75) 



2 



dQ 

However, a thermal bath does carry energy. Hence the energy that is attributed to a jet in a thermal medium will 
depend on the scale used to probe the jet. 

In order to analyze the effect of evolution on the fragmentation functions at finite temperature, we must first 



resolve the difficulty mentioned in Section III A: It is not truly consistent to treat the fragmentation of a low- 
energy parton without taking the fragmentation of the bath itself into account. To separate the contribution from 
distinct partons to our evolution equations from that of the bath, we introduce the 'assimilation functions' Pq,g{z/T) 
and pq_g{z/T) = 1 — pq_g{z /T). Phenomenologically, pq(^g^{z /T) represents the probability that a quark (gluon) of 
momentum z will ther malize with a bath of temperature T before it has time to fragment. Keeping track of only the 



first term in Eq. (74a), we write 

Q^^Fy,/q{y; X, Q2) dy/x = Q2_1_ ^ (y, dy);Q^) 



Stt Jq X 
a,(Q2) dz 



asiQ ) r ^^^^(^,.^)p^^,/T)P{z ^ {y,dy);Q') 

o , %q{z;x}pq{z/T)P{z ^ {y,dy);Q'^) . (76) 

ZTT Jq X 

Here, P{x — > {y,dy);Q^) is the probability that a quark of momentum x will decay into a hadron with momentum 
between y and y + dy and transverse momentum k'j^ < Q^. In the first term on the right-hand side of the second 
equality, the quark maintains its identity after the bremsstrahlung. Hence the use of F/^^q^y; z,Q'^) to describe the 
subsequent fragmentation is appropriate. In the second term, however, the quark thermalizes with the medium before 
the fragmentation takes place. Therefore, it is actually the thermal bath that produces the final-state hadron. 

Although the latter contribution naturally appears when calculating the radiative corrections to the fragmentation 
functions at finite temperature as defined in Eq.(p2|), it is inconsistent with the idea of a distinct parton fragmenting 
into an observed hadron while propagating through a thermal bath. Experimentally, it belongs in the background 
and should be removed before the data are analyzed. For this reason, we retain only the first term in our evolution 
equation. This procedure does not come without a price. In particular, the fragmentation function described by our 
evolution equation is no longer rigorously associated with the thermal matrix element appearing in (|2^). However, 
this truncation is necessary in order to describe a phenomenologically meaningful function. Discrepancies between the 
truncated fragmentation function described below and that defined in Eq.(|2^) are important only when we consider 
hadrons which are easily generated by the bath. Since these hadrons are difficult to distinguish experimentally from 
the background, their contribution does not concern us here. 

Extending the above analysis to the full evolution equation dTJ), we arrive at 



Q'^ -^Fh/q{y;x,Q'^) = -^-^p- "7 [lqq{z]x)pq{z /T)Fh/q{y; z,Q'^) ^ -iqg{z]x)pg(z /T)Fh/g{y] z,Q'^)\ (77a) 

0^ '^^h/giv, X, Q^) = ^'^^ ^ y \lgq{z; x)pq{z /T) F^/ siv, z, Q^) + lgg{z] x)pg{z /T)Fh/ g{y] z, Q^)] . (77b) 

As explained above, these evolution equations contain contributions only from processes in which the parton propa- 
gating through the thermal medium remains independent of the medium. 
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The assimilation functions, Pq^g{z/T), are easy to understand phenomenologically, but difficult to define rigorously. 
They can be determined from phenomenological models of parton thermalization, such as those used in Ref . fl^ . All 
that is required for Eq. ( |77| ) to be well-defined mathematically is that Pq{C) faster than ( and Pg{C) ~^ faster 
than as C ^ 0, along with the obvious restriction < Pg,g(C) < 1- In addition, we impose the physical restriction 
Pq^giO ^ as C ^ oo. 

At this point, we are ready to investigate the effects of evolution on simple observables of the parton jet. Defining 
the jet energies 



Eqix, Q^)^J2 Tl^ F^/q{y; x, Q^) ; 
, Jo 2; X 



we have 



d 



zdz 



^^^Eg{x,Q') 



dQ^ 



27r Jq X 
c^siQ^J f°° zdz 

X 



27r 







2 \lqq{z\ x)pq{z/T)Eq{z, Q^) + -fqg{z] x)pg{z /T)Eg{z, Q^)] ; 

,2 [lgq(.z;x)pqiz/T)Esiz,Q^) +^gg{z;x)pg{z/T)Eg{z,Q^)] . 



(78a) 
(78b) 

(79a) 
(79b) 



These functions represent the fraction of the initial parton energy carried by the final hadron jet. Given initial 
conditions, we can solve this set of equations numerically to obtain the scale dependence of the jet energy. 

We can get some idea of the effect of evolution at large x ^ T by considering Eq^g{z, Qq) = 1 at some scale Qq. 
This is precisely the zero-temperature value of £'g,g, so it should be approximately valid at high energy. This ansatz 
leads to 



dQ- 



■Eq{x,Q^) 



dQ' 



where 



Eg{x,Q^) 



2tt 



f^^^{x)+tf^{x) 
2nfm{x)+f^^Kx) 



— (-) lpp'{z]x)pp,{z/T) 
X \x/ 



(80a) 
(80b) 

(81) 



Contributions to Eq.( pOaD can be understood in terms of the energy loss and gain to different partonic degrees of 
freedom due to the basic bremsstrahlung/absorption process. The first and second terms represent energy gain to the 
quark and gluon components of the fragmenting quark, respectively. Since the vacuum does not carry any energy, 
energy losses to the quark degrees of freedom must be compensated completely by gains in the gluon sector at zero 
temperature. Explicit calculation yields 



(82) 



in agreement with our expectations. 

As the temperature is increased from zero, several effects work to alter this result. To analyze their contributions 
separately, it is useful to replace ^^ ^(z/r) with 1 — pq g{z/T) in rqq,g{x). The '1' contains the vacuum contributions 
described above, as well as the thermal effects of induced emission and absorption to Tqq (x) and those of Pauli 
blocking and quark annihilation to Tqg (x) . The assimilation functions take the energy flow from explicit parton 
degrees of freedom to the heat bath into account. These contributions act to reduce the amount of energy available 
to both quark and gluon degrees of freedom in the jet. 

Induced gluon emission and absorption allow the thermal bath to change explicitly the energy of the distinct parton 

— (2) 

degrees of freedom. These processes only involve thermal gluons, so they do not contribute to Tqg (x). The form of 
the splitting function Pqq{zf) guarantees that the absorption process wins out over that of induced emission, giving 
rise to the net fractional energy gain 
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{{T/xf) 



(83) 



to the quark degrees of freedom. The analogous processes of Pauh blocking and quark annihilation affecting the 
gluonic degrees of freedom cancel almost exactly, leaving a net energy gain from annihilation that is suppressed by 
exp(— x/r). 

Our distinct parton degrees of freedom lose energy to the medium mainly through the direct bremsstrahlung process. 
In order to obtain an explicit expression for this contribution to the energy loss, we must take a specific form for 
the assimilation functions. For illustrative purposes, we choose the hard cut-off Pq.g{C) — (^iVq,g ~ 0) where rjq^g ~ 1 
governs the thcrmalization of partons introduced to our thermal medium. This choice yields the fractional energy 
losses 



-0{{T/x) 

3 



52f(^J{x) = -2Cp ^ + Cf -Icf + OieM% ^ ^/T)) 



(84a) 
(84b) 



The appearance of a linear (in T/x) term in Eq. (84b) is quite striking since it dominates all other corrections as 
X oo. This term is generated by the divergent behavior of the soft gluon emission rate. Since Pqg{zf) — > as 
— > 0, the fractional energy carried by soft gluons saturates. This causes the amount of energy absorbed by the 
bath to be proportional to the size of the energy window that thermalizes with the bath, which is dictated by the 
temperature T . Such a linear term will appear in the asymptotic expansion (as x ~f oo) for any choice of assimilation 
function satisfying the requirements as we have given before. 

Putting all of this together and performing similar calculations in the gluon sector, we rewrite Eq.(pQ) as 



dQ 



dQ' 



■Eg{x,Q^ 



27r X 



2Vg 



2x / 



VqT 

2x 



l7r'- + Oi{T/xr) 

D X 



(85a) 



2n X 



2Ca \Vg -Vg 



VgT 



2nfTF Vq 



VqT 

2x 



I X 



— TT — 

3 X 

2 



OiiT/xY) 



(85b) 



As before, the linear term in the gluon sector is generated by the divergence of the soft gluon bremsstrahlung rate. 
These expressions indicate that the parton energy which contributes to a hadronic jet is reduced by an amount 
approximately proportional to the temperature of the thermal medium as the transverse momentum cut-off of the jet 
is increased. The effect is slightly smaller at lower energies due to absorption processes. This phenomenon arises from 
the fact that increasing the momentum cut-off increases the amount of transverse phase space available to partonic 
sub-processes, which increases the amount of energy lost to the bath. Although our result was derived from a constant 
input distribution and rather simplistic assimilation functions, the qualitative effect of evolution is expected to be 
similar to that of a more realistic calculation. In particular, the dominant behavior of the thermalization mechanism 
over the phenomena of induced emission, absorption, annihilation and Pauli blocking is expected to persist in a 
complete calculation. 

Fluctuations in the thermal medium can also screen the quark flavor of a jet. The valence structure of a jet is 
measured by the quark number 



Nq{x,Q^) 



E 

Hq) 



dy 



[Fh/q{y; X, Q^) - Fh/q{y; x, Q^)] 



(86) 



where the sum goes over all hadrons h whose valence structure includes q. Only fragmentation processes in which 
the leading quark contributes to the valence structure of the observed hadron can contribute to Eq. (p6|), since all 
others will cancel in the difference. In vacuum, Nq{x, Q^) represents the net flavor of the jet. Since QCD interactions 
conserve flavor, 



Nq{x,Q^)\T^^ = \ 



(87) 



is independent of scale or parton energy. At finite temperature, the leading quark in a jet can scatter into the thermal 
bath or annihilate with an anti-quark in the bath. As before, these processes generate a dependence of the quark 
number on both scale and parton energy. Eq.(77a) leads to the evolution equation 
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dQ 



27r 



dz . 



■ 799(2^; x)pq{z/T)Nq{z, Q^) 



(88) 



for Nq{x,Q''). 

As before, we illustrate the effect of this evolution equation at high energy by considering Nq{x,QQ) = 1 at the 



scale Qo- This gives 



dQ- 



■Nqix,Q^) 



2n 



r(V(x) 



Using the same assimilation function as above, we arrive at the expression 



ni\x) = -CF 



T 



2x 4 x 



0{{T/xf) 



(89) 



(90) 



which indicates a net transfer of flavor from the jet to the bach as the momentum scale is increased. It is interesting 
to note that while the effect of the medium on energy loss becomes larger with the energy of the leading parton, 
its effect on flavor loss diminishes as the energy is increased. This is essentially because the fractional rather than 
absolute energy loss is comparable to the net flavor loss. 

We can also consider the energy loss of a quark propagating through a thermal medium from a slightly different 
point of view. Our interpretation of Eq. (fz^ ) indicates that 



dPqq{x z, 



dz dQ^ 



2n 7.,(^;x)p-,(z/T)- — 



(91) 



represents the probability that a leading quark of energy x will split into a quark with energy between z and {z + dz) 
and transverse momentum between and + dQ^ and a gluon. The effective energy loss of the leading quark can 
be summarized as the expectation value of the emitted gluon's energy. 

To obtain a number for this quantity, we must first determine a range of values for the probing scale. Since we 
have identified phenomenologically with the transverse momentum k'^ between the quark and gluon, it is natural 
to use the kinematic limit 



fci < 4z|x- 



(92) 



as shown in Ref. This limit depends on the details of the process, so our result cannot be directly related 
to a renormalization group analysis. For this reason, the appearance of the running coupling in Eq. (|9l| ) is not fully 
consistent with the idea of renormalization. We will take the coupling evaluated at some fixed scale, Q^Tcharacteristic 
of our process. The infrared scale /z^, given as the Debye screening mass, is taken as a lower limit of the fc^ integral. 
These considerations lead to the quark fractional energy loss 



as{Ql)CF 



dz 



kl 



lqq{z;x)pq{z/T) 



2ti 



4 , 4x2 49 
-log — - — 
3 /i^ 18 



^ flog-^'^''''^ 



log 



4xT 
1^ 



2 - 7B + 




0{{Tlxf) 



(93) 



where we have again used the hard cut-off form of Pq{z/T). 

The first term in this expression represents quark energy loss due to vacuum bremsstrahlung, and is independent 
of the medium. Some of the energy lost to the quark through these processes re-appears in the gluon fragmentation, 
while the rest is lost to the medium. The difference is unimportant here since we do not track the gluon degrees of 
freedom. The second term represents the shift in net quark number of the jet as the scale varies. Its sign indicates that 
it represents energy gain to the quark when the energy is large. This is because we have considered the normalized 
energy loss to one quark, while the net quark number of the jet is actually less than one due to annihilation effects in 
the medium. The result is that the energy per quark increases. The third term corresponds to energy that comes from 
the medium directly either from absorption/emission processes or via soft quark thermalization. This contribution 
indicates that the absorptive processes in the plasma overcome the emissive processes in the absence of thermalization 
{r]q = 0), leading to a net energy gain. The details of the calculation depend on our choice of assimilation function, 
as well as transverse momentum cut-offs, but the results are nonetheless interesting. We note here that secondary 
scattering effects, such as those studied in Ref. appear only at the next-to-leading-log approximation. These are 
not taken into account in our formulae. 
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V. CONCLUSION 



In this paper we have taken a first step in the study of parton fragmentation functions in a thermal QCD medium via 
thermal field theory. We extended the definition of the parton fragmentation functions to the case of a thermal medium 
with finite temperature which will reduce to the conventional fragmentation functions in the zero temperature limit. 
We then calculated the radiative corrections to these fragmentation functions within the framework of perturbative 
QCD at finite temperature and derived the corresponding DGLAP evolution equations. The modified evolution 
equations can be written in the same form as at zero temperature. However, the effective (or medium modified) 
splitting functions depend explicitly on the temperature and contain terms corresponding to various physical processes 
such as radiation, absorption, backward and forward scattering, etc. The introduction of the temperature causes the 
effective splitting functions to depend explicitly on the value of the initial parton energy. This complicates the 
renormalization group analysis of the effective fragmentation functions and causes us to modifiy their intepretation 
slightly. However, we have shown that the modified evolution equations are self-consistent. There are both linear and 
logarithmic infrared divergences in the radiative corrections. Some cancel among themselves while others are canceled 
by the virtual corrections. 

As an example of the application of the evolution equations, we also derived the evolution equations for the net 
quark number and the energy loss (gain). We find that one has to introduce 'assimilation functions' in order to 
separate the hadronization of the thermal bath from that of the initial parton jet. We find similarly as in an earlier 
study Q that the absorption of thermal gluons wins out over the stimulated emission, due to the form of q qg 
splitting function. This effective medium-induced energy gain will then reduce the total energy loss experienced by a 
propagating parton. 

As we have emphasized in the Introduction, we have considered only the first order corrections in the leading 
logarithmic approximation. In this case we only included sequential parton emission or absorption. Collision-induced 
bremsstrahlung occurs beyond this order of perturbation theory. The intermediate thermal gluon will no longer be 
on- shell and its propagator will also develop an imaginary part. This will result in collision-induced radiation. When 
the formation time of the gluon radiation is larger than the coUisional rate which is contained in the imaginary part 
of the propagator, there will be LPM suppression of the induced radiation We hope to address this problem in 
future studies. 

Acknowledgements. — This work was supported by the the Director, Ofhce of Energy Research, Office of High 
Energy and Nuclear Physics, Division of Nuclear Physics, and by the Office of Basic Energy Science, Division of 
Nuclear Science, of the U.S. Department of Energy under Contract No. DE-AC03-76SF00098, and by National 
Natural Science Foundation of China under projects 19928511, 10135030 and 10075031. E. W. thanks LBNL Nuclear 
Theory Group for its hospitality during the completion of this work. X.-N. W. thanks A. Majunder for helpful 
discussions. E. W. and X.-N. W. also thank the Department of Physics, Shandong University, where part of the work 
was completed, for its hospitality 



[1] 

[2] 

[3] 
[4] 
[5] 

[6] 

[7] 

[8] 

B 



H. Mueller (Ed.), Perturbative Quantum Chromodynamics, World 
438 (1972); Yu. L. Dokshitzer, Sov. Phys. JETP 46, 641 (1977); G. 



J. C. Collins, D. E. Soper, and G. Sterman, in: A 
Scientific, Singapore, 1989, and the references therein. 
V. N. Gribov and L. N. Lipatov, Sov. J. Nucl. Phys. 15 
Altarelli and G. Parisi, Nucl. Phys. B126, 298 (1977). 
S. Bethke J. Phys. G26, R27 (2000) | |arXiv:hep-ex/000402l{ . 

E. Wang and X.-N. Wang, Phys. Rev. Lett. 87, 142301 (2001) |arXiv:nucl-th/0106043|. 

X.-N. Wang, Z. Huang and I. S arcevic, Phys. Rev. Let t. 77, 231 (1996) |arXiv:hcp-ph/9605213[ ; X. N. Wang and Z. Huang 



Phys. Rev. C 55, 3047 (1997) [ [a.rXiv:hep-ph/9701227[ . 
M. Gyulassy and X.-N. Wang, Nucl. P hys. B 420, 583 (1994) ]arXiv:nucl-th/9306003tl ; X.-N. Wang, M. Gyulassy and M. 



Pliimer, Phys. Rev. D 51, 3436 (1995) [ |irXiv:hep-ph/9408344[ . 

R. Baier, Y. L. Dokshitzer, A. H. Mueller, S. Peigne a nd D. Schiff, Nucl. Phys. B483, 291 (1997) ]arXiv:hep-ph/9607355 1 



Nucl. Phys. B 484, 265 (1997) ]arXiv:hep-ph/9608322[ 
M. Gyulassy, P. Levai a nd I. Vitev, Nucl. Phys. B594, 371 (2001) | |arXiv:nucl-th/0006010t ; Phys. Rev. Lett. 85, 5535 (2000) 
I arXiv:nucl-th/0005032| . 

U. A. Wiedemann, Nucl. Phys. B 588, 303 (2000) ||arXiv:hep-ph/000512S|]. Nucl. Phys. A 690, 731 (2001) |arXiv:hep- 



ph/0008241| 



[10] X.-N. Wang, Phys. Lett. B 485, 157 (2000) |arXiv:nucI-th/0003033| ] . 

[11] L. D. Landau and I. J. Pomeranchuk, Dolk. Akad. Nauk. SSSR 92, 92(1953); 



A. B. Migdal, Phys. Rev. 103, 1811 (1956). 



20 



[12] 

[13] 
[14] 
[15] 
[16] 
[17] 
[18] 
[19] 



A. H. Mueller, Phys. Rev. D18, 3705 (1978); S. Gupta and A. H. Mueller, Phys. Rev. D20, 118 (1979); A. H. Mueller, Phys 
Rep. 73, 237 (1981). 

R. D. Field, in: Application of Perturbative QCD, Frontier in Physics Lecture, Vol. 77, Addison Wesley, 1989. 
D. McLerran and T. Toimela, Phys. Rev. D31, 545 (1985). 
C. Collins, D. E. Soper and G. Sterman, Nucl. Phys. B261, 104 (1985). 
C.Collins and J. Qiu, Phys. Rev. D39, 1398 (1989). 
A. Weldon, Phys. Rev. D28, 2007 (1983). 
Geiger and B. MuUer, Nucl. Phys. B 369, 600 (1992). 

Arnold, G . D. Moore and L. G. Yaffe, JHEP 0111, 057 (2001) |arXiv:hep-ph/0109064[ ; JHEP 0206, 030 (2002) |arXiv:hep- 



l/0204343|. 



21 



